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PREFACE 


Part  of  the  RAND  research  program  consists  of  basic 
supporting  studies  in  mathematics.  The  mathematical 
research  in  this  Memorandum  is  concerned  with  techniques 
for  obtaining  approximate  solutions  to  systems  of  nonlinear 
differential  equations,  which  arise  in  a  variety  of  physical 
problems . 
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SUMMARY 

In  some  recent  papers,  [l],  [2],  [3],  we  have  applied' 
various  closure  techniques  to  the  problem  of  obtaining 
approximate  solutions  of  nonlinear  differential  equations. 
These  questions  fall  within  the  area  of  differential 
approximation.  Starting  with  a  vector  system 

dx  ,  v 
dt =  g(x)’ 

we  replace  it  by  another  system 
«  =  h(y)’ 

with  simpler  analytic  and  computation  properties  in  such 
a  way  that  II  x  -  y  il  is  small.  In  some  cases  we  are 
content  merely  to  have  h(y)  linear,  i.e.,  of  the  form 
Ay  +  b;  in  other  cases,  we  want  the  dimension  of  y  to 
be  considerably  less  than  that  of  x. 

In  this  paper  we  wish  to  consider  the  case  where  the 
original  system  is  infinite  dimensional  and  the  approximating 
system  is  to  be  finite  dimensional.  Infinite  dimensional 
systems  arise  in  a  natural  fashion  from  the  consideration 
of  partial  differential  equations. 
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CLOSURE  TECHNIQUES  FOR  INFINITE  SYSTEMS 
OF  DIFFERENTIAL  EQUATIONS 

1 .  INTRODUCTION 

In  some  recent  papers,  [l],  [2],  [3],  we  have  applied 
various  closure  techniques  to  the  problem  of  obtaining 
approximate  solutions  of  nonlinear  differential  equations. 
These  questions  fall  within  the  area  of  differential 
approximation.  Starting  with  a  vector  system 

d-D  ^  =  g(x), 

we  replace  it  by  another  system 

(1.2)  ^  =  h(y), 

with  simpler  analytic  and  computation  properties,  in  such 

a  way  that  II  x  -  y  IJ  is  small.  In  some  cases  we  are 

content  merely  to  have  h(y)  linear,  i.e.,  of  the  form 
Ay  +  b;  in  other  cases,  we  want  the  dimension  of  y  to 
be  considerably  less  than  that  of  x. 

In  this  paper  we  wish  to  consider  the  case  where  the 
original  system  is  infinite  dimensional  and  the  approximating 
system  is  to  be  finite  dimensional.  Infinite  dimensional 
systems  arise  in  a  natural  fashion  from  the  consideration 
of  partial  differential  equations.  For  example,  if  we 
look  for  a  solution  of 


(1.3) 
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in  the  form 

(1.4)  u(x,  t)  =  Z  u.  (t)  eikx , 

k  K 

we  obtain  an  infinite  system  of  ordinary  differential 
equations  for  the  components  u^(t);  see  [4]. 

The  problem  of  closure  is  now  more  difficult  since 
there  is  no  immediate  way  of  expressing  the  components 
corresponding  to  higher  harmonies  in  terms  of  the  components 
for  smaller  k.  In  this  paper  we  shall  discuss  two  approaches. 

2.  EXTRAPOLATION  AND  SELF-CONSISTENCY 

Let  us  suppose  that  we  have  an  infinite  linear  system 
of  the  form 

du.  oo 

(2J)  dt~  =  ^  ak{  u{  -  \<°>  *  ck  -  k  =  l,  2,  ... 

For  a  discussion  of  existence  and  uniqueness  of  solutions 
to  equations  of  this  type,  and  an  examination  of  when  the 
sections 

duk  N 

(2'2)  dtT  =  ^  aki  uf  >  uk(o)  =  ck,  k  =  1,  2,  . . N, 

yield  solutions  which  converge  to  the  solutions  of  (2.1), 
see  [5] . 

Instead  of  using  the  closure  technique  implicit  in 
(2.2),  namely,  uk  =  0,  k  >  N  +  1,  we  wish  to  approximate 
to  the  remainder  by  means  of  linear  combinations  of  the 


preceding  uk, 
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N 

(2.3)  Z  1.1  =  E  b,  t  u,  , 

l  >  N  +  1  k{  t=l  kl  1 


see  [l],  [2],  [3], 

The  coefficients  bk^  are  to  be  chosen  so  that 


r  N  2 

(2-4)  J  [  £  a.  f  ii-  -  I  b  f  uJ  dt 

0  i  >  N+l  £  {=1  k*  * 


is  a  minimum.  We  obtain  in  this  way  the  linear  algebraic 
equations  for  the  b^ 


(2.5) 


l  > 


l  +1\i  l  "I  Ur  «  ■  £  brf  /  Ur  « 


r  =  1,  2,  .  .  . ,  N. 


The  usual  difficulty  now  confronts  us.  How  do  we  compute 


the  integrals  J  u-  u  dt  involving  the  unknown  solution? 
0  *  r 

Consider  first  the  case  where  1  <  r,  f  <  N.  To 

obtain  these  integrals,  we  use  the  finite  system  in  (2.2). 

Call  the  solutions  u^°\  The  coefficients  of  b  *  in 

r  rf 


(2.5)  are  then 


(o) 


dt .  Observe  that  these 


quantities  can  be  computed  directly  in  the  course  of 

obtaining  the  u^0^  by  adjoining  to  (2.2)  the  equations 

♦ 

dw 


(2.6) 


r  i 


dt 


=  ur  u|  >  wrf(o)  =  0  ’ 


and  asking  only  for  the  values  wf^(T). 
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The  more  difficult  problem  is  that  of  the  calculation 

of  /.  u.dt  for  r=l,  2,  w,  £  >  N  +  1 .  Here 

0  r  £  - 

we  use  extrapolation  techniques.  Keep  r  fixed  and  le,t  £ 

vary  over  the  integers,  £  =  1,  2,  ...  It  is  reasonable 

to  suspect  that  f ,  =  J  u  u,  dt  will  be  a  well-behaved 

£,  r  g  r  £ 

sequence.  Consequently,  if  we  possess  the  values 

f 1  r>  *2  r ’  • • • »  r’  we  can  use  an7  a  number  °f 
extrapolation  techniques  (see  [6])  to  obtain  the  values  of 

f *  „  for  £  >  N  +  1 . 

The  system  in  (2.5)  then  has  the  form 


(2.7) 


E 

M  >  £  >  N  +  1 


ak£  f£,r 


N 

E 

£=1 


(T)  . 


Here  M  is  a  cut-off  number,  such  as  2N,  which  depends 
upon  the  size  of  the  coefficients  a^  and  the  rapidity 
of  convergence  of  the  infinite  series. 

Solving  (2.7)  numerically,  we  obtain  the  coefficients 
br£^  ’  We  use  the  suPerscri-Pt  to  indicate  the  fact  that 
these  are  the  first  approximations.  To  obtain  higher 
approximations,  we  use  self-consistency  techniques. 

In  place  of  (2.2),  let  us  now  use 


(2.8) 


dt 


N 

E 

£=1 


u, 


+ 


N 

E 

£=1 


Call  the  solutions  of  this  equation  u 
We  now  proceed  as  before  to  calculate 


u£  ’ 

(1) 

k  ’ 


uk<0)  =  ck. 


k  =  1, 
(T),  f 


(1) 

£,  r 


N. 


and 
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bkl  '■  With  the  new  coefficients  b^ ,  we  introduce  the 
equation 

du,  N  N  /•,  \ 

(2'9)  ar  ’  £  ak!  ui  +  ^  bki  V  uk(0)  =  V 

This  process  is  repeated  until  the  values  of  the  b,  t 
settle  down. 

3.  NONLINEAR  EQUATIONS 

In  the  case  where  the  infinite  system  arises  from  a 
nonlinear  equation  such  as  (1.3),  we  can  obtain  a  finite 
system  in  a  simple  fashion.  Write 

(3.1)  u  =  E  u  (t)elkx  , 

N  iki  <  N  k 

replace  u  u^  by  u^Cu^)^,  and  in  this  fashion  obtain  the 
higher  harmonics  in  terms  of  the  lower  harmonics . 
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